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I.  INTRODUCTION 


The  purpose  of  this  investigation  is  to  analyze  the  protec¬ 
tion  afforded  by  active  Beacon  Collision  Avoidance  System  (BCAS) 
alert  logics  in  view  of  certain  anticipated  compromises  that 
must  be  made  in  the  logic  parameters.  It  has  been  argued  that 
such  compromises  are  necessary  in  order  to  prevent  the  alert 
rate  that  the  system  would  generate  in  congested  traffic  areas 
from  becoming  too  high. 

An  attempt  has  also  been  made  in  this  investigation  to  use 
real  air  traffic  data  to  test  the  assumptions  on  which  the  anal¬ 
ysis  is  based.  However,  scope  limitations  have  made  it  neces¬ 
sary  to  rely  upon  data  that  already  exist  in  accessible  form. 

The  source  chosen  for  this  purpose  is  a  briefing  presented  in 
February  1980  (Ref.  1)  by  the  Mitre  Corporation,  which  is  cur¬ 
rently  responsible  for  developing  active  BCAS  alert  logics  for 
the  FAA.  Reasons  for  choosing  this  particular  data  source  are 
discussed  in  Chapter  III. 

Not  surprisingly,  in  the  data,  which  are  derived  from 
Automated  Radar  Terminal  System  (ARTS)  tapes  containing  65  hours 
of  radar  recording  at  Houston  International  Airport,  there  is 
no  instance  of  a  mid-air  collision.*  Because  such  events  are 
indeed  rare,  it  Is  not  feasible  to  use  the  incidence  of  actual 
collisions  as  a  measure  of  BCAS  protection  in  a  test  involving 
real  air  traffic.  However,  if  it  Is  conceded  that  some  minimum 
separation  between  aircraft  is  necessary  for  safe  flight,  then 
a  minimum  miss  distance  standard  can  be  imposed  and  its 

^According  to  Ref.  2,  fatalities  due  to  mid-air  collisions  are 
only  0.2  percent  of  all  passenger  fatalities  due  to  transpor¬ 
tation  accidents  in  general. 


maintenance  adopted  as  the  basis  for  empirically  assessing  the 
amount  of  protection  that  a  collision  avoidance  system  can  fur¬ 
nish.  This  principle,  the  use  of  which  was  already  implied  in 
Ref.  1,  in  fact  defines  what  is  meant  here  by  the  protection 
that  a  given  version  of  BCAS  alert  logic  will  provide. 

A  tool  that  is  fundamental  in  achieving  the  stated  goals 
of  this  investigation  is  a  mathematical  model  of  the  statistical 
distributions  of  air  traffic  density  and  velocity  over  the 
Houston  terminal  region.  A  similar  analytical  characterization 
was  originally  devised  (Ref.  3)  to  fit  the  same  statistical 
distributions  in  the  Los  Angeles  Standard  Traffic  Model  (Ref.  4), 
which  consists  of  data  synthetically  generated  to  simulate 
traffic  conditions  that  were  predicted  for  the  Los  Angeles 
Basin  in  1982.  Chapter  III  and  Appendix  A  furnish  details  con¬ 
cerning  the  consistency  of  this  traffic  model  with  various 
empirical  results  taken  from  Ref.  1. 

The  evaluation  of  active  BCAS  protection  occurs  in  Chapter 
IV.  Among  the  questions  considered  there  are:  (1)  basic  fac¬ 
tors  that  should  be  taken  into  account  in  the  selection  of 
active  BCAS  alert  logic  parameter  values;  (2)  for  a  given  choice 
of  parameters.  If  no  alert  occurs,  the  probability  that  a  pre¬ 
scribed  miss  distance  standard  will  be  violated;  (3)  if  an  alert 
does  occur,  the  probability  that  it  is  unnecessary,  i.e.,  that 
the  miss  distance  standard  would  not  be  violated  if  the  alert 
were  ignored;  (4)  If  an  alert  occurs,  the  probability  distribu¬ 
tion  of  the  time  available  for  the  aircraft  to  react  and  per¬ 
form  an  avoidance  maneuver;  (5)  loss  of  protection  due  to  the 
use  of  slant  range  in  the  horizontal  alert  algorithm  instead  of 
horizontal  range. 

Implicit  in  the  discussion  of  protection  in  Chapter  IV  is 
the  likelihood  that  an  aircraft  pilot’s  perception  of  the  qual¬ 
ity  of  alerts  generated  by  the  BCAS  logic  may  be  as  important 
as  the  actual  quantity.  That  is,  it  seems  reasonable  to  suppose 


that  no  one  would  object  to  frequent  alerts  if  they  appear  to 
be  justified  by  circumstances;  however,  if  the  overwhelming 
majority  of  alerts  are  viewed  as  unnecessary,  a  pilot  may  reject 
the  system  even  when  the  alert  rate  is  moderate.  It  also  seems 
reasonable  to  expect  an  even  greater  loss  of  confidence  in  the 
system  if  it  should  produce  many  unnecessary  alerts  and  also, 
on  occasion,  fail  to  produce  one  that  would  have  been  perceived 
as  justified. 

The  most  important  conclusions  resulting  from  the  present 
investigation  are  summarized  in  Chapter  V.  Although  a  number 
of  such  conclusions  were  reached,  the  primary  objective  has 
been  to  formulate  analytical  procedures  that  may  be  useful  in 
selecting  active  BCAS  logic  parameters.  For  example,  with  the 
aid  of  such  procedures  it  may  be  possible  in  the  future  to 
incorporate  into  the  system  adaptive  BCAS  parameter  selection 
(desensitization),  tailored  both  to  the  observed  traffic  dis¬ 
tribution  and  each  individual  aircraft  according  to  its  aero¬ 
dynamic  or  other  relevant  characteristics. 
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II.  DESCRIPTION  OF  ACTIVE  BCAS* 


A.  GENERAL  CHARACTERISTICS 

Active  BCAS  (Beacon  Collision  Avoidance  System)  is  a  coop¬ 
erative  airborne  collision  avoidance  system  of  the  beacon-trans¬ 
ponder  typ'e.  It  will,  in  principle,  protect  a  user  against  any 
other  aircraft  equipped  with  an  Air  Traffic  Control  Radar  Beacon 
System  (ATCRBS)  transponder  and  altitude  encoder.  It  will  also 
protect  against  future  participants  in  the  Discrete  Address 
Beacon  System  (DABS)  that  is  intended  to  replace  ATCRBS. 

In  its  operation  BCAS  relies  upon  measurements  of  an  intrud¬ 
er's  relative  range  and  altitude.  It  also  requires  the  intruder* 
range  and  altitude  rates,  which  it  obtains  by  tracking. 

BCAS  generates  the  measured  data  by  transmitting  the  same 
ATCRBS  interrogation  signals  that  are  usually  transmitted  by 
ground  installations.  An  intruder  will  respond  in  the  usual 
manner  with  altitude  (and  identity)  encoded  reply  signals.  The 
signal's  delay  time  provides  the  range  and  its  encoded  message 
the  altitude. 

This  data  is  processed  to  determine  whether  a  threat  exists. 
Then  the  BCAS  conflict  resolution  logic  determines  what  commands 
or  advisories  are  needed  for  protection  of  the  aircraft. 


*The  Information  contained  in  this  chapter  was  extracted  from 
Ref.  5  and  Ref.  6  and  reinterpreted  slightly  to  emphasize  a 
point  of  view  that  is  appropriate  to  the  present  investigation. 
The  alphanumeric  variable  and  parameter  names  used  in  this 
and  subsequent  chapters  are  those  given  in  Table  2-1,  Ref.  5 • 


B.  THREAT  CRITERIA 


BCAS  employs  two  types  of  threat  criteria:  (1)  horizontal, 
which  depends  on  the  intruder's  relative  range  and  range  rate; 
(2)  vertical,  which  depends  on  the  intruder's  relative  altitude 
and  altitude  rate.  A  threat  is  recognized,  and  an  alert  is 
generated,  only  if  both  criteria  are  satisfied  simultaneously 
and  if,  by  means  of  a  linear  projection  of  the  altitude  rate, 
it  is  estimated  that  the  vertical  separation  of  the  aircraft 
will  be  less  than  a  prescribed  threshold,  called  ALIM,  when  the 
aircraft  reach  their  position  of  closest  approach. 

The  horizontal  threat  criterion  depends  upon  two  fixed 
parameters,  RDTHR  and  HI,  and  two  prescribed,  but  flexible, 
parameters:  TRTHR,  designated  here  by  r,  and  DMOD,  designated 
here  by  R  .  The  vertical  threat  criterion  depends  upon  a  single 
prescribed  flexible  parameter,  TVTHR,  which,  however,  is  always 
assigned  the  same  value  as  TRTHR,  and  a  fixed  parameter  ZTHR. 

For  a  horizontal  separation  R  and  separation  rate  R,  if 
the  aircraft  are  closing  at  a  rate  that  is  greater  than  the 
fixed  parameter  RDTHR,  the  horizontal  threat  criterion  is  the 
condition 

R  +  t  R  £  R_  . 


-RDTHR  £  -R  £  RDTHR 

R  is  replaced  by  -RDTHR.  Otherwise,  for  aircraft  that  are 
separating,  the  criterion  becomes 

R  R  <  HI  . 

The  vertical  criterion  is  satisfied  if  either  of  two  con¬ 
ditions  are  met.  For  a  relative  altitude  A  and  altitude  rate  A, 
the  first  condition  is 


A  +  x  A  n  0, 


and  the  second  is 


£ 


A  *  ZTHR. 


C.  DESENS ITIZATION 

3ecause  of  traffic  conditions,  if  nothing  were  done  to 
alter  the  threat  logic,  it  is  expected  that  an  aircraft  enter¬ 
ing  terminal  air  space  would  experience  increasingly  larger 
BCAS  alert  rates.  In  addition,  the  close  proximity  of  aircraft 
in  normal  terminal  flight  patterns  would  tend  to  cause  unneces¬ 
sary  and  disruptive  alerts. 

Therefore,  a  provision  has  been  made  in  the  system  design 
concept  to  make  the  threat  logic  more  restrictive  by  altering 
the  flexible  threat  criterion  parameters  when  BCAS  operates 
in  certain  predetermined  zones.  This  process,  known  as  desensi¬ 
tization,  will  be  implemented  automatically  by  means  of  radar 
signal  transmissions  from  certain  strategically  located  ground- 
based  radar  beacon  transponders  (RBXs).  An  RBX  may  also  pro¬ 
vide  certain  useful  information,  such  as  the  local  altitude 
relative  to  sea  level. 

At  the  time  that  this  report  was  prepared,  five  different 
levels  of  desensitization  could  be  imposed.  The  following  rules 
were  then  current  for  aircraft  below  10,000  ft.  In  level  1  the 
BCAS  alert  logic  is  inhibited,  as  is  its  interrogation  of  RBXs. 
In  level  2  the  alert  logic  is  inhibited  but  RBX  interrogation 
is  permitted.  In  level  3  the  value  of  DMOD  ( RQ )  is  set  at  0.1 
nmi,  t  (TRTKR  and  TVTHR)  is  set  at  18  sec  when  the  intruder  is 
equipped  with  BCAS  and  20  sec  when  it  is  not,  and  the  value  of 
ALIM  is  set  at  3^0  ft.  In  level  4  the  value  of  DMOD  is  set  at 
0.3  nmi,  x  is  set  at  25  sec,  and  the  value  of  ALIM  is  set  at 
3^0  ft.  In  level  5  the  value  of  DMOD  is  set  at  1.0  nmi,  t  is 
set  at  30  sec,  and  the  value  of  ALIM  is  set  at  440  ft. 

All  aircraft  above  10,000  ft  are  in  level  5*  For  such 
aircraft,  however,  altitude  thresholds  with  larger  values  re¬ 
place  ALIM  because  of  increased  altimeter  error. 
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III.  AIR  TRAFFIC  STATISTICS 


A.  DATA 

From  the  outset  of  this  investigation  it  was  felt  that  some 
empirical  data  would  be  needed  to  support  the  conclusions 
reached  by  analysis  that  otherwise  must  depend  largely  upon  a 
set  of  unsubstantiated  assumptions.  Unfortunately,  most  of  the 
presently  available  air  traffic  data  are,  at  least  in  part, 
synthetic  in  nature. 

For  example,  NAFEC#  has  performed  a  series  of  experiments 
(Refs.  7-9)  to  determine  alert  rates  and  the  effect  of  BCAS  on 
pilots  and  air  traffic  control  operations.  These  experiments 
have  involved  real  people;  however,  all  are  based  on  computer 
generated  traffic  data  designed  to  simulate  conditions  that 
are  thought  to  be  specific  to  certain  airports. 

Other  potentially  useful  sources  of  data  also  have  an 
artificial  character.  Among  these  is  the  source  that  undoubtedly 
furnishes  the  most  detailed  air  traffic  data  that  is  presently 
available — the  Los  Angeles  Standard  Traffic  Model  (Ref.  4). 

None  of  the  data  that  comprise  that  model  were  actually  derived 
from  direct  observation  of  real  air  traffic. 

There  are  some  references  (Refs.  1,  10,  11)  that  do,  at 
least,  summarize  certain  real  air  traffic  statistics.  Of  these. 
Ref.  1  contains  the  most  extensive  statistical  information. 

For  this  reason  this  report  relies  upon  Ref.  1  for  its  empirical 
data  base.  Various  statistical  summaries  from  that  source  are 
used  here  to  test  the  foregoing  theoretical  analysis  as  a  means 
of  verifying  indirectly  the  assumptions  on  which  the  analysis 
depends . 

^National  Aviation  Facilities  Experimental  Center. 


Reference  1  contains  processed  and  reduced  data  extracted 
from  ARTS  tapes  covering  a  65-hour  period  of  radar  observation 
at  Houston  Intercontinental  Airport.-  The  processing,  by  com¬ 
puter,  involved  the  application  of  BCAS  alert  logic,  for 
various  choices  of  algorithm  parameters,  to  the  Houston  traf¬ 
fic  in  order  to  determine  how  the  choice  of  parameters  might 
affect  the  alert  rate  and  the  protection  that  the  logic  affords 

Of  particular  interest  for  present  purposes  are  the 
results  contained  in  Table  1,  which  gives  the  number  of  alerts 
that  occurred  for  several  combinations  of  tau  and  DMOD  values 
for  two  different  altitude  threshold  standards  (ALIM).  Table  1 
consists  of  the  relevant  part  of  a  similar  but  more  extensive 
table  appearing  in  Ref.  1. 

TABLE  1.  DISTRIBUTION  OF  ALERTS  (PARTIAL  DATA) 


DMOD 


0.3  n mi 
0.3  nmi 
0.1  nmi 
0.1  nmi 
0.1  nmi 
0.1  nmi 


ALIM 


370  ft 
340  ft 
370  ft 
340  ft 
370  ft 
340  ft 


Positive  Commands 


25  sec 
25  sec 
20  sec 
20  sec 
1 8  sec 
18  sec 


Reference  1  also  contains  some  useful  statistical  data  in 
the  form  of  histograms  showing  the  distribution  of  speeds  and 
separations  of  aircraft  in  conflict.  In  addition,  L.  Zarrelli 
kindly  supplied  some  histograms  that  did  not  appear  in  the  orig 
Inal  presentation.  One  of  these  is  also  relevant  to  the  pur¬ 
poses  of  this  report:  a  histogram  giving  the  observed  distri¬ 
bution  of  closing  speeds  of  aircraft  in  conflict.  These  histo¬ 
grams  are  reproduced  in  Appendix  A  (Figs.  A-l,  A-3  and  A-4). 
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B.  THE  THEORETICAL  TRAFFIC  MODEL 


Reference  3  indicates  how  the  distribution  of  low-altitude 
air  traffic  density  about  a  point  of  reference  in  the  Los 
Angeles  Basin  and  the  distribution  of  aircraft  velocities, 
defined  by  the  Los  Angeles  Standard  Traffic  Model  data  (Ref.  4), 
can  be  represented  analytically.  The  mathematical  expressions 
used  for  this  purpose,  although  simple,  fit  the  data  with  an 
accuracy  that  is  consistent  with  the  statistical  fluctuations 
of  the  several  traffic  samples  provided  in  Ref.  4. 

This  mathematical  model  of  air  traffic  can  be  formulated 
as  follows.  The  horizontal  coordinates  of  an  aircraft  picked 
at  random,  located  by  means  of  a  cartesian  coordinate  system 
with  its  origin  at  the  traffic  center,  are  independent,  iden¬ 
tically  distributed  Gaussian  random  variables.*  Each  coordinate 
has  a  zero  mean  and  the  same  standard  deviation  aR.  Similarly, 
the  cartesian  components  of  the  aircraft  velocity  vector  are 
also  independent,  identically  distributed  Gaussian  random  vari¬ 
ables  with  zero  mean  and  the  same  standard  deviation  aQ. 

One  consequence  of  these  assumptions  is  that  the  aircraft 
density  is  Rayleigh  distributed  in  range  R  from  the  traffic 
center.  That  is,  the  probability  density  PR  of  aircraft  dis¬ 
tribution  over  range  is  given  by 


f Actually,  it  was  found  that  the  distribution  is  bimodal. 
However,  since  the  fraction  of  aircraft  in  one  mode  is  small 
compared  to  that  in  the  other,  the  single  mode  approximation 
described  here  is  still  quite  accurate. 


Another  consequence  is  that  the  separation  R  between  a  pair  of 
aircraft  selected  at  random  is  also  Rayleigh  distributed  with 

SR  =  aR/T  » 

replacing  crR  in  the  definition  of  P^.  A  third  consequence  is 
that  the  components  of  the  relative  velocity  between  the  two 
aircraft  are  also  independent,  identically  distributed  Gaus¬ 
sian  random  variables  with  mean  zero  but  with  a  standard  devia¬ 
tion 


°o  /T 


From  this  it  follows  that  the  closing  speed  is  Rayleigh  distri¬ 


buted,  i.e.,  with  a  probability  density  Pv  given  by 


P  = 


2o 


It  was  determined  in  Ref.  3  that  a  good  fit  Is  obtained 
for  the  Los  Angeles  Standard  Traffic  Model  by  setting 

-  20  nmi 


and 


122  ft/sec  (72.2  kn)  . 


Since  no  detailed  information  is  available  for  the  distribution 
of  traffic  and  aircraft  velocities  at  Houston  it  is  not  possible 
to  verify  directly  a  similar  fit  of  the  model  to  the  observed 
Houston  traffic.  However,  the  statistical  summaries  provided 
by  Ref.  1  are  available,  at  least,  to  test  the  compatibility 
of  such  a  model  with  the  underlying  data  recorded  on  the  ARTS 


C.  APPLICATION  OF  THE  THEORETICAL  TRAFFIC  MODEL  TO  THE  HOUSTON 
DATA 

In  order  to  apply  the  theoretical  traffic  model  to  the 
Houston  data  given  in  Ref.  1,  it  is  necessary  to  calculate  the 
appropriate  alert  probabilities  in  terms  of  the  assumed  values 
of  x  and  Rq  (DMOD).  This  is  done  in  Appendix  A. 

There  it  is  shown  that  the  four  quantities  t,  R  ,  a„,  and 

O  K 

ay  can  be  subsumed  into  two  independent  dimensionless  param¬ 
eters  p  and  k,  defined  by 


v 

3T' 


An  inspection  of  these  equations  reveals  that  (1)  p  is  inde¬ 
pendent  of  the  aircraft  density;  (2)  k  is  independent  of  DMOD; 
(3)  a  change  in  r  is  equivalent  to  a  change  in  the  aircraft 
velocity  distribution;  (4)  the  product  pic  is  independent  of  the 
aircraft  velocity  distribution  and  t. 

Figure  A- 2  illustrates  the  result  of  one  type  of  calcula¬ 
tion  done  in  Appendix  A.  It  provides  curves  of  the  probability 
P(k,p)  that  two  co-altitude  aircraft  chosen  at  random  wij.1  be 
in  conflict,  i.e.,  will  satisfy  the  horizontal  alert  condition 
(Eq  1  in  Chapter  IV,  Section  A),  as  a  function  of  p  for  various 
values  of  k.  It  Is  assumed  that  the  alert  rate,  or  the  number 
of  alerts  occurring  during  the  period  covered  by  the  Ref.  1 
ARTS  tapes,  is  proportional  to  this  probability  for  the  appro¬ 
priate  values  of  p  and  k. 

In  order  to  fit  the  theory  to  the  histograms  from  Ref.  1 
it  is  assumed  that  virtually  all  conflicts  occur  at  altitudes 
below  10,000  ft  where  the  theoretical  traffic  model  is  appli¬ 
cable.  It  is  also  implicit  In  the  analysis  that  the  histo¬ 
grams  are  based  on  a  random  sample  from  the  traffic.  This  is 
only  approximately  true  for  the  histogram  depicting  the  distri¬ 
bution  of  speeds  of  aircraft  in  conflict,  since  the  fact  that 


each  such  aircraft  is  one  of  a  pair  for  which  the  alert  con¬ 
dition  is  satisfied  tends  to  bias  the  sample  against  low-speed 
aircraft . 


Figures  A-l,  A-3  and  A-4  illustrate  the  results  of  other 
calculations  done  in  Appendix  A,  to  determine  the  theoretical 
distribution  that  best  fits  each  histogram.  The  calculations 
indicate  that  the  aQ  value  of  72.2  kn  appropriate  to  the  Los 
Angeles  Standard  Traffic  model  is  too  small.  The  figures  con¬ 
tain  curves  based  on  crQ  values  of  83-3  kn,  which  would  cor¬ 
respond  to  99  percent  confidence  that  no  aircraft  would  violate 
the  250  kn  speed  limit,  and  106  kn,  which  appears  to  provide 
the  best  fit.  The  theoretical  distributions  are  essentially 
independent  of  the  aircraft  density  variation,  since  varying  tc 
from  0.01  to  0.1  produces  only  insignificant  modifications  in 
the  curves. 

Using  the  probability  P(<,p)  illustrated  by  the  curves  of 
Fig.  A-2,  Appendix  A  tests  the  ability  of  the  theory  to  predict 
the  alert  changes  due  to  changes  in  t  and  DMOD  indicated  in 
Table  1.  This  is  done  for  the  three  possible  aircraft  velocity 
distributions  characterized  by  oQ  values  of  72.2  kn,  83.3  kn, 
and  106  kn.  The  results  are  given  in  Table  2. 

TABLE  2.  ALERT  REDUCTION  DUE  TO  CHANGES  IN  x  AND  DMOD 


Change  in 

x  From  25  Actual  %  Predicted  % 

sec  to  ALIM  Change  in  Change  in  Alerts _ Error  (%) 


(sec) 

(ft) 

Alerts 

ITT 

83.3 

106 

72.2 

83.3 

i 

■ 

20 

370 

54.1 

55.1 

53.2 

49.9 

1.0 

0.9 

4.2 

20 

3  40 

53.2 

55.1 

53.2 

49.9 

1.9 

0.0 

3.3 

18 

370 

59.7 

62.4 

61.2 

58.7 

2.7 

1.5 

1.0 

18 

60.1 

62.4 

61.2 

58.7 

2.3 

1.1 

1.4 
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According  to  Table  2  the  best  predictions  occur  for  oQ  = 
83.3  kn.  However,  the  prediction  error  is  small  for  all  three 
aQ  values,  including  106  kn,  which  provides  the  best  theoretical 
fit  to  the  histogram  data.  Thus,  it  appears  that  predictions 
of  the  change  in  the  alert  rate  due  to  changes  in  the  values 
of  t  and  DMCD  is  insensitive  to  the  aircraft  velocity  distribu¬ 
tion. 

Also,  as  pointed  out  in  Appendix  A,  the  change  is  virtually 
independent  of  the  aircraft  density  distribution.  This  obser¬ 
vation  leads  to  the  natural  conjecture  that  changes  in  t  and 
DMOD  would  have  the  same  effect  on  fractional  changes  in  the 
alert  rate  at  any  other  airport . 
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IV.  ALERT  LOGIC  PROTECTION 


A.  PARAMETER  SELECTION 

For  the  purposes  of  the  foregoing  analysis  it  is  assumed 
that  collision  alerts  are  triggered  by  the  idealized  horizontal 
algorithm  embodied  in  the  condition 

R  +  x  R  *  R0,  (1) 

where  R  is  the  separation  between  encountering  aircraft,  R  is 
the  separation  rate,  and  t  and  RQ  (termed  DMOD  in  BCAS  parlance), 
are  prescribed  positive  parameters. 

As  shown  in  Ref.  12,  the  algorithm  (1)  will  only  furnish 
protection  against  collisions  under  all  circumstances  when  the 
DMOD  parameter  RQ  is  at  least  large  enough  to  satisfy  the  rela¬ 
tion 


m 


i 


R0  *  +  |  OT2  . 

In  this  relation  is  a  bound  on  the  horizontal  separation 
measurement  error,  is  a  bound  on  the  horizontal  separation 
rate  measurement  error,  U  is  a  bound  on  the  maximum  anticipated 
relative  acceleration  of  the  encountering  aircraft,  and  t  is 
the  effective  time  required  to  respond  to  an  alert  and  maneuver 
to  a  safe  vertical  separation.  If  a  guaranteed  minimum  miss 
distance  requirement  is  to  be  imposed  as  well,  the  miss  distance 
standard  D  must  also  be  added  to  the  right  side  of  the  inequality. 

Allowance  for  protection  against  acceleration  is  sometimes 
discounted  because  it  is  usually  assumed  that  aircraft  only 
accelerate  during  turning  maneuvers  and  that  only  ten  percent  of 
aircraft  in  a  terminal  region  are  expected  to  be  turning  at  any 


given  time.  However,  for  some  aircraft  an  allowance  for  accel¬ 
eration  may  always  be  necessary. 

In  performing  a  vertical  avoidance  maneuver  after  an  alert, 
acceleration  must  be  applied  vertically,  but  it  may  not  be  pos¬ 
sible  for  some  aircraft  to  create  an  acceleration  vector  that 
does  not  also  have  a  horizontal  component.  When  this  is  the 
case,  that  component  must  be  taken  into  account  in  the  hori¬ 
zontal  collision  algorithm  like  any  other  horizontal  accelera¬ 
tion. 

The  danger  of  ignoring  this  effect  is  similar  to  that  in 
the  situation  discussed  in  Ref.  14,  wherein  a  fast  aircraft  is 
overtaking  a  slower  one  in  a  tail  chase  and  the  avoidance  maneu¬ 
ver  is  supposed  to  be  horizontal.  There  it  was  assumed  that 
acceleration  takes  place  normal  to  the  flight  path  throughout 
the  maneuver,  so  that  the  aircraft  maintains  constant  speed  and 
turns  with  a  constant  angular  rate.  It  was  found  that,  because 
forward  velocity  along  the  original  flight  path  decreases,  if 
the  slower  aircraft  is  the  one  actempting  to  perform  the  avoid¬ 
ance  maneuver  a  collision  may  be  unavoidable  no  matter  how  large 
the  estimated  time  to  collision  may  be. 

A  parameter  similar  to  DMOD  should  be  defined  for  the 
vertical  collision  algorithm,  for  which  BCAS  presently  employs 
an  unmodified  tau  criterion,  since,  at  the  very  least,  altimeter 
measurement  errors,  particularly  the  error  in  estimating  alti¬ 
tude  rates,  are  far  from  negligible.  Moreover,  when  there  is 
a  non-zero  altitude  rate  the  existence  of  a  vertical  accelera¬ 
tion  is,  if  anything,  more  likely  than  the  existence  of  a 
horizontal  acceleration  in  general. 

In  the  selection  of  collision  algorithm  parameters  it 
seems  logical  to  choose  t  first  because,  unless  RQ  is  at  least 
as  large  as  the  minimum  allowable  miss  distance  D,  the  protection 
afforded  by  DMOD  depends  upon  r.  On  the  other  hand,  if  the  value 
of  t  is  reduced  without  changing  the  DMOD  value  the  protection 
afforded  by  DMOD  will  appear  to  be  greater. 

18 


However,  this  will  be  true  only  for  those  aircraft  that 
can  respond  to  an  alert  quickly  enough  to  make  use  of  the  infor¬ 
mation  that  it  provides.  On  the  other  hand,  reducing  DMOD 
without  changing  t  reduces  the  protection  for  all  aircraft, 
although  it  is  also  true  that  the  protection  afforded  by  any 
DMOD  value  is  greater  for  the  aircraft  that  can  respond  more 
quickly  to  an  alert. 


Some  compromise  in  protection  will  undoubtedly  occur  even 
for  level  5*,  for  which  a  DMOD  value  of  1  nmi  and  a  t  value  of 
30  sec  have  been  selected.  For  example,  even  if  no  allowance 
is  made  for  measurement  error  or  a  minimum  miss  distance,  the 

p 

acceleration  bound  U  for  which  the  term  1/2  Ut  contributing  to 
Rq  is  equal  to  1  nmi  is  just  13*5  ft/sec2  or  about  0.42  g.  If 
a  not  unreasonable  allowance  of  1,000  ft  is  made  for  the  total 
measurement  error  contribution  and  1,000  ft  for  the  miss  dist¬ 
ance,  the  system  will  protect  against  accelerations  up  to  about 
0.3  g.  For  aircraft  that  can  respond  quickly  enough  to  accom¬ 
modate  a  t  value  of  25  sec,  the  protection  against  acceleration 
without  taking  into  account  measurement  error  or  a  minimum  miss 
distance  is  about  0.6  g.  With  a  total  measurement  error  plus 
miss  distance  allowance  of  2,000  ft,  the  protection  would  still 
take  into  account  accelerations  up  to  0.4  g.** 

The  parameter  t  should  be  a  sum  of  three  terms:  (1)  the 
data  updating  interval  which,  because  of  BCAS  logic  rules  (Ref. 
5),  If  an  alert  is  not  missed,  will  be  2  or  3  secs;  (2)  the 
pilot-aircraft  response  time  which,  according  to  the  results  of 

^Reference  4  defines  BCAS  protection  levels  1-5  that  result  from 
desensitization  of  the  system  in  air  terminal  regions,  rang¬ 
ing  from  no  CAS  operation  at  level  1  to  maximum  protection  at 
level  5  (see  p.  7) . 

**In  fact,  if  1,000  ft  is  an  appropriate  allowance  for  measure¬ 
ment  error  when  t  is  30  sec,  the  measurement  error  term  would 
lessen  when  t  is  25  sec;  therefore,  the  protection  is  actually 
somewhat  greater  than  this  estimate  indicates. 


a  NAFEC  experiment  (Ref.  15)  is  not  likely  to  exceed  8  sec, 
and  (3)  the  time  t  required  for  a  vertical  avoidance  maneuver. 

In  estimating  t^  it  will  be  assumed  that  the  avoidance 

e  2 
maneuver  begins  with  a  constant  vertical  acceleration  (8  ft/sec  , 

according  to  Ref.  5)  which  is  applied  until  a  limiting  altitude 

rate  v^  is  reached.  Reference  5  prescribes  16.67  ft/sec  for 

v^  in  a  climbing  maneuver  and  25  ft/sec  in  a  descending  maneuver. 

The  assumption  of  a  constant  vertical  acceleration  may  be 
incorrect.  For  example,  if  the  motion  is  like  that  of  a  hori¬ 
zontal  turn  in  which  the  angular  turning  rate  is  constant, 
although  the  magnitude  of  the  acceleration  may  have  a  constant 
value  a,  its  vertical  component  would  be  a  function  of  time, 

a 

namely,  a  cos  — t,  where  v  is  the  aircraft's  speed  when  it 
begins  the  maneuver.  The  estimate  for  tg  should  be  modified 
appropriately  if  the  vertical  acceleration  is  not  constant, 
e.g.,  by  assuming  that  a  is  an  equivalent  average  value  for  an 
acceleration  that  may  be  considerably  larger  at  its  peak. 

The  choice  of  t  must  be  at  least  consistent  with  the 
e 

case  of  a  co-altitude  aircraft  encounter  with  no  relative  alti¬ 
tude  rate.  In  the  case  of  non  co-altitude  encounters  for  which 
the  aircraft  are  closing  vertically  the  choice  of  tQ  depends 
upon  the  maximum  altitude  rate  expected  for  an  intruder.  As 
this  expected  rate  for  an  intruder  approaches  the  appropriate 
avoidance  maneuver  limit,  16.67  ft/sec  or  25  ft/sec,  specified 
for  the  BCAS  logic,  the  required  value  of  t  becomes  arbitrarily 
large  unless  both  aircraft  perform  cooperative  maneuvers. 

For  the  co-altitude  case  a  vertical  acceleration  a  occurs 
for  a  time  tQ  given  by 


i.e.,  acceleration  occurs  for  2.1  sec  in  a  climbing  maneuver 
and  3.1  sec  in  a  descending  maneuver.  Since  the  value  of  t 
will  surely  be  greater  than  either  of  these  times,  the  limiting 


speed  of  16.67  ft/sec  or  25  ft/sec  will  always  be  reached 
before  the  maneuver  is  completed.  The  vertical  separation  z 
when  the  maneuver  is  completed  is  given  by 


z  *  V1  fce  -  IT  ■ 

When  the  level  5  value  of  t,  i.e.,  30  sec,  is  used  t  is 
19  sec.  Then  calculations  of  z  show  that  a  vertical  separation 
of  299  ft  will  be  achieved  in  a  climbing  maneuver  and  ^36  ft  in 
a  descending  maneuver.  For  a  t  of  25  sec  t  is  14  sec  and  the 
separations  will  be  1 66  ft  and  236  ft  respectively. 

For  a  final  separation  of  just  160  ft,  in  a  climbing  maneu¬ 
ver  it  would  be  unnecessary  to  use  a  value  greater  than  10  sec 
for  t  .  This  would  correspond  to  a  r  value  of  only  21  sec. 

B.  ALGORITHM  SENSITIVITY  TO  PARAMETER  CHANGES 

The  quantity  Rq  is  essentially  an  allowance  for  measure¬ 
ment  error  and  any  relative  acceleration  due,  in  general,  to 
possible  turning  maneuvers  by  the  encountering  aircraft.*  The 
quantity  t  is  an  allowance  for  the  amount  of  time  that  would 
be  required  for  an  aircraft  to  perform  a  vertical  avoidance 
maneuver  after  an  alert.* 

If  these  interpretations  of  RQ  and  t  actually  govern  their 
choice  their  values  will  be  such  as  to  guarantee  a  successful 
avoidance  maneuver  in  any  encounter  for  which  the  assumed  bounds 
on  acceleration  and  escape  maneuver  time  are  satisfied.  However, 
if  for  some  values  of  either  parameter  the  condition  (1)  is  not 
adequate,  it  is  also  possible  to  revise  the  algorithm  to  ensure 
that  a  safe  avoidance  maneuver  can  be  accomplished  by  increasing 
the  value  of  the  other  parameter  sufficiently.  That  is,  if  t  is 
insufficient  by  an  amount  At,  for  example,  for  a  particular  R  the 
corrective  effect  would  be  the  same  if,  instead  of  increasing  t 


»Cf.  Ref.  12  and  Ref.  13- 
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by  At,  Rq  were  changed  to  RQ  -  At  R.  This  would  amount  to  an 

increase  of  -AtR  in  DMOD  (or,  conversely,  the  change  in  t  would 
AR 

be  -  — when  the  aircraft  are  closing.  Therefore,  the  equiva- 
R 

lent  change  in  DMOD  would  be  large  for  large  closing  speeds, 
small  for  small  closing  speeds,  and  it  would  actually  be  neces¬ 
sary  to  decrease  the  value  of  DMOD  in  order  to  obtain  the  equiv¬ 
alent  effect  when  the  encountering  aircraft  are  separating. 

This  suggests  that  it  is  dangerous  to  draw  conclusions  about 
the  effectiveness  of  a  particular  choice  of  RQ  and  t  when  the 
evidence  is  based  on  a  statistical  sample  of  events,  whether 
simulated  or  real,  unless, extreme  caution  is  exercised  in  the 
interpretation  of  the  events. 

C.  THE  RISK  WHEN  THERE  IS  NO  ALERT 

As  observed  in  Ref.  4,  since  midair  collisions  are  quite 
rare  it  will  be  easier  in  practice  to  evaluate  a  collision 
avoidance  system  in  terms  of  the  near  misses,  that  is,  viola¬ 
tions  of  prescribed  separation  standards,  that  the  CAS  prevents 
rather  than  actual  collisions  that  it  might  have  prevented. 
Although  the  usual  interpretation  of  (1)  (Refs.  12  and  13)  refers 
to  the  avoidance  of  a  collision  rather  than  a  violation  of  a 
horizontal  separation  standard,  it  is  possible,  nevertheless, 
by  making  RQ  sufficiently  large,  to  include  in  the  algorithm  a 
limit  on  the  minimum  miss  distance  permitted  in  an  encounter. 

In  order  to  examine  the  relation  between  the  collision 
algorithm  associated  with  (1)  and  the  minimum  miss  distance 
standard.  It  is  convenient  to  use  a  coordinate  system  of  the 
type  introduced  by  J.  Holt  in  Ref.  12  for  describing  the  hazard 
region  defined  by  (1).  The  origin  of  the  coordinate  system  is 
located  at  one  of  two  encountering  aircraft  with  a  coordinate 
axis  parallel  to  the  relative  velocity  vector  of  the  other. 

The  other  aircraft  is  located  by  polar  coordinates  relative  to 
this  axis,  which  is  assumed  to  be  oriented  in  a  direction  oppo¬ 
site  that  of  the  relative  velocity.  Assuming  a  constant 
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relative  velocity  of  magnitude  V,  the  separation  rate  is  given 
in  terms  of  the  angular  coordinate  0  by 

R  -  -V  cos  6.  (2) 

A  cross  range  velocity  of  magnitude  VQ  is  also  implied,  given 

by 

V0  »  V  sin  6.  (3) 

The  condition  (1)  is  satisfied,  then,  at  all  points  within  the 
cardioid  region  bounded  by  the  curve 

R  *  Rq  +  tV  cos  9,  (4) 


illustrated  in  Fig.  1. 

In  the  analyses  that  follow,  for  the  most  part,  it  is  to 

be  understood  that  R  has  had  removed  from  it  the  measurement 

o 

error  allowance,  and  either  linear  flight  is  assumed  or  the 

acceleration  allowance  has  also  been  removed.  The  error  allow' 

ance  R„„  has  the  form 
oe 


where  is  a  bound  on  the  error  in  measuring  R  and  is  a 

bound  on  the  error  in  measuring  R.  The  acceleration  allowance 

R„„  has  the  form 
oa 


Roa  "  I  U'[2> 

where  U  is  a  bound  on  the  magnitude  of  the  possible  relative 
acceleration. 


In  addition  to  polar  coordinates  it  is  also  convenient  to 
use  cartesian  coordinates  defined,  as  usual,  by 

x  =  R  cos  8,  y  =  R  sin  0. 

For  linear  flight,  by  definition,  the  intruding  aircraft  located 
at  (x,y)  moves  In  the  direction  of  decreasing  x  and  so  that  y 
remains  constant.  Thus,  the  minimum  miss  distance  is  always 
equal  to  y. 


23 


In  discussing  levels  of  safety  afforded  by  the  condition  (1) 
perhaps  the  most  basic  question  that  should  be  addressed  is 
what  separation  is  assured  if  no  alert  at  all  occurs  during  an 
encounter.  In  terms  of  the  cardioid  illustrated  in  Pig.  1,  this 
can  happen  only  when  the  ordinate  y  associated  with  the  intru¬ 
der’s  position  is  larger  than  the  maximum  ordinate  ym  of  points 
on  the  curve  determined  by  (M). 

The  maximum  ordinate  ym  is  a  solution  of  the  equation 


dR 


dy  _  dff 


sin  8  +  R  cos  8 


dx  cos  0  -  R  sin  6 


0, 


which  is,  more  explicitly,  equivalent  to 

-tV  sin20  +  (Rq  +  tV  cos  6)cos  0 

The  equation  (5)  is  equivalent  to 

cos2  0  +  R^  cos  0  -  i  =  0 
o  d 

2tV 


for  which  possible  solutions  are 

-  R 


=  0. 


(5) 


cos  0  = 


±  ^Rq2+8t2V2 
4tV 


(6) 


Then,  since 


ym  =  R  sin  0  *  (Rq+tVcos  0)sin  0, 


on  substituting  from  (6)  it  can  be  seen  that  ym  is  given  by  the 
greater  of 


r±  -  5  P  V  V 


R02  +  8T2V2)(i  + 


’o s  4 


R  ^+8t2V2> 
o 
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> 
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It  can  be  seen  from  (7)  that  y  Is  an  increasing  function 

in 

of  V  as  V  varies  from  0  to  Thus,  while  ym  can  be  arbitrarily 

large,  the  smallest  value  that  it  can  have  is  Rq.  Therefore, 
if  there  is  no  alert  a  minimum  miss  distance  of  at  least  RQ  is 
guaranteed. 

The  minimum  miss  distance  ym  given  by  (7)  as  a  function  of 

velocity  V  can  be  illustrated  by  a  single  curve,  shown  in  Fig.  2, 

if  y  is  measured  in  units  of  R^  ft  and  V  is  measured  in  units 
■'m  o 

of  Rq/t  ft/sec.  It  can  be  seen  by  inspection  of  Fig.  2  that 
the  curve  is  essentially  linear  for  all  values  of  V  above  5 
and  deviates  only  slightly  from  linear  for  values  between  0 
and  5. 

The  curve  illustrates  the  fact  that  even  when  R  is  quite 

o 

small,  say  5  to  10  ft,  although  It  is  very  possible  that  a  normal 
separation  standard  such  as  1,000  ft  might  be  violated,  it  is  / 
still  likely,  at  least,  that  the  miss  distance  will  exceed  100 
ft.  For  example.  If  Rq  is  5  ft  and  t  is  25  sec  a  V  of  40  cor¬ 
responds  to  a  relative  velocity  of  8  ft/sec.  According  to  Fig.  2 
the  value  of  ym  for  V  equal  to  40  is  about  20.7,  which  is  equiv¬ 
alent  to  more  than  103  ft.  In  other  words,  as  long  as  the 
relative  velocity  is  at  least  8  ft/sec  the  miss  distance  will 
be  greater  than  103  ft.  For  a  x  of  20  sec  the  velocity  corre¬ 
sponding  to  40  would  be  10  ft/sec;  since  this  is  greater  than 
8  ft/sec,  a  miss  distance  exceeding  103  ft  would  be  slightly 
less  likely. 

A  probability  distribution  for  the  minimum  miss  distance 


would  give  a  more  precise  Idea  of  the  risk  associated  with  a 
collision  algorithm  when  an  alert  condition  does  not  occur  in 


wmmmm 


an  encounter.  This  probability  distribution  depends  upon  the 
knowledge  of  a  probability  distribution  for  the  relative  veloc¬ 
ity,  which  may  be  a  property  of  the  local  air  traffic. 


For  the  Los  Angeles  Standard  Traffic  model  (Ref.  4),  in 
the  lower  altitudes,  each  velocity  component  of  an  aircraft 
appears  to  be  normally  distributed  with  a  mean  of  zero  and  a 
standard  deviation  of  about  122  ft/sec.  This  would  imply  a 
similar  distribution  for  each  component  of  the  relative  veloc¬ 
ity  vector  between  two  aircraft,  but  with  a  standard  deviation 
of  122  2  »  173  ft/sec.  The  relative  velocity  magnitude  V 

would  then  be  Rayleigh  distributed  with  a  standard  deviation  of 
173  ft/sec.  Thus, 


where , 


P(V  *  VQ) 


in  this  case. 


VdV  = 


1 


(8) 


a  =  173  ft/sec. 


Equation  (8)  used  in  conjunction  with  (7)  provides  the  proba- 
bility  distribution  P(ym  *  D)  of  interest. 

The  complimentary  probability 

P(ym  >  D>  =  1  -  p(yra  *  D) 

provides  a  lower  bound  on  the  probability  that  the  miss  distance 

will  be  at  least  D  if  there  is  no  alert.  Figure  3  contains 

curves  of  this  bound  for  normalized  values  o  of  a (measured 

in  units  of  R  A  ft/sec)  and  D  (measured  in  units  of  ft) 
o  o 

covering  the  CAS  avoidance  algorithm  parameters  that  are  likely 
to  occur  at  altitudes  below  10,000  ft.  The  value 


a  =  2.37 

for  one  curve  shown  in  Fig.  3  corresponds  to  a  a of  173  ft/sec 
when  RQ  is  0.3  nmi  and  t  is  25  sec. 


If  the  Rq  value  of  0.3  nmi  is  assumed  to  cover  measurement 
errors  for  which  a  total  allowance  of  1,000  ft  is  to  be  made, 
the  curve  labelled  with  the  normalized  value 

a  =  5.25 

provides  the  appropriate  probability  bound.  The  measurement 
error  allowance  of  1,000  ft  might,  for  example,  result  from  a 
range  error  of  150  ft  and  a  range  rate  error  of  slightly  more 
than  20  kn.  An  examination  of  the  curve,  which  would  then  cor¬ 
respond  to  a  reduced  value  for  Rq  of  824  ft,  shows  that  the 
probability  of  at  least  a  1,000  ft  miss  distance  is  98.8  percent 
or  better. 

D.  THE  RISK  WHEN  AN  ALERT  OCCURS 


Except  for  the  case  of  pop-up  targets,  an  alert  occurs 
when  an  intruder  is  on  or  near  the  boundary  of  the  cardioid 
shown  in  Fig.  1.  A  reasonable  measure  of  safety  provided  by  the 
CAS  algorithm  (1)  might  then  be  the  amount  of  time  tD  available 
for  maneuver  after  such  an  alert  before  the  separation  R  between 
the  encountering  aircraft  is  reduced  to  some  miss  distance  D 


chosen  as  a  minimum  separation  standard. 

In  calculating  tn  it  may  be  observed  that,  because  of  (2), 


y-y-vv 


so  that  when  an  alert  occurs,  implying  at  that  instant  the  rela¬ 
tion 


R  +  t  R  =  R 


(1*) 


in  accordance  with  (1), 


When  the  aircraft  separation  is  D 


I 


-*D=V^v. 


(11) 


since  it  is  assumed  that  the  flight  is  linear  and,  therefore, 
that  y  is  constant.  It  follows  from  (10)  and  (11)  that 


4 


1 

D2  -  R2 

(R-R0)2' 

u  n 

TJ  V2 

(12) 


If  xA  is  the  horizontal  distance  of  the  intruder  -to  the 
point  of  closest  approach  at  the  time  when  the  alert  occurs, 
then  the  time  tA  between  the  alert  and  the  closest  approach  is 
given  by 


.  fA  ,  R  cos  6 

ZA  V  V 


Then,  because  of  (!’)  and  (2),  it  follows 


that 


R ( R-R  ) 

4.  _  _ o 


tV 


The  definitions  of  tA  and  x^  imply  that 


(13) 


(14) 


V 

£ 


t  as  +■  _ 

rD  ZA  V  • 

Thus,  because  of  (14)  and  (12)  it  follows  that 


R(R-R  )  \D2-1 

t  =  _ 2 _ " 


/d2  r2 

L  Hn 

1  D  -n 

— i 

CM 

> 

CM 

H 

±_l 

-  tV*"  V 

It  may  happen,  however,  that  when  the  alert  occurs 


(15) 


in  which  case  has  no  meaning  since  the  separation  will  always 
be  at  least  D.  Because  of  (9)  the  condition  (16)  is  equivalent 


/r2-d‘ 


When  the  aircraft  are  closing,  which  is  the  only  non-trivial 
case,  this  condition  is 


V  2  -RR 


2  2 
R  -D* 


R(R-R0) 
tJ  r2-d2 


(17) 


because  of  (1* ) . 


The  relation  (17)  defines  V  .  which  is  the  largest  value 
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that  V  can  have  when  there  is  any  risk  that  the  separation 
standard  can  be  violated  at  all.  The  minimum  value  V .  ol 


of  V 


is,  of  course,  the  magnitude  -R  of  the  separation  rate;  i.e.. 


(18) 


According  to  (17)  and  (18),  then. 


rj — 5  vmin 

Jr  -d* 


(19) 


From  (15)  it  follows  that  the  corresponding  times  tDmin  and  tDma 
to  violation  of  the  separation  standard  are  given  by 


'Dmin 


'Dmax 


'Drain 


(“) 


(20) 


(21) 


The  quantity  tDmin  given  by  (20)  is  the  least  possible 
amount  of  time  available  for  maneuvering  before  the  separation 
standard  is  violated  after  an  alert.  It  corresponds  to  the  case 


of  a  head-on  collision.  The  quantity  tpmax  given  by  (21)  is 
the  greatest  amount  of  time  available  for  maneuvering  when  the 
alert  is  necessary,  i.e.,  if  the  separation  standard  will,  in 
fact,  be  violated  unless  an  avoidance  maneuver  is  performed. 

It  corresponds  to  the  case  in  which  the  miss  distance  will  be 
exactly  equal  to  the  separation  standard.  Figures  4a  and  b 
illustrate  the  functional  dependence  of  tDmin  and  tDmax  on  the 
separation  R  for  various  values  of  D. 


The  solid  line  curves  in  Fig.  4a  correspond  to  the  param 

eters  of  the  previous  example  in  which  RQ  is  0.3  nmi  (with  a 

1,000  ft  measurement  error),  t  is  25  sec  and  D  is  1,000  ft. 

Since  in  Fig.  4a  tD  is  normalized  to  units  of  t  sec  and  R  and 

D  to  units  of  R„  ft,  for  this  case  the  normalized  value  of  D 
o 

is  1.2136  as  indicated. 


For  R  very  large  compared  to  Rq  and  D,  (20)  and  (21)  show 
that  the  time  available  to  perform  a  maneuver  to  ensure  the 
required  separation  is  t,  which  is  the  same  as  the  time  provided 
by  the  CAS  algorithm  to  avoid  an  actual  collision.  This,  of 
course,  is  not  a  particularly  surprising  result. 

If  Rq  Is  larger  than  D,  it  has  already  been  seen  that  the 
separation  requirement  is  guaranteed  when  no  alert  occurs. 

On  the  other  hand.  If  an  alert  does  occur,  according  to  (20),  as 
the  range  R  approaches  RQ  the  time  available  to  perform  a  maneu¬ 
ver  that  guarantees  the  required  separation  increases,  and,  in 
fact,  becomes  arbitrarily  large.*  If  Rq  is  less  than  D  it 
follows  from  (20)  and  (21)  that  the  time  available  to  maneuver 
before  the  separation  standard  Is  violated  is  always  less  than 
t  unless  the  separation  is  sufficiently  large,  i.e.,  unless 


R  >  D 
R 


(22) 


Such  an  occurrence  may  still  not  produce  absolute  safety  (cf. 
Sec.  A). 


MAX 


The  probability  distribution  of  t^  after  an  alert  has 
occurred  can  be  used  to  assess  the  risk  of  a  CAS  algorithm. 
This  probability  distribution  is  directly  related  to  the  dis¬ 
tribution  of  the  relative  cross-range  velocity  component  VQ, 
which,  consistent  with  the  assumptions  stated  earlier,  will  be 
taken  to  be  a  normally  distributed  random  variable  with  mean 
zero. 

The  standard  deviation  a  for  V0  is  f~2  times  that  asso¬ 
ciated  with  a  velocity  component  of  a  single  aircraft.  Hence, 
for  example,  for  V0  in  the  Los  Angeles  Standard  Traffic  Model 
o  would  be  173  ft/sec. 

Prom  (2)  and  (3)  it  follows  that 


2  *2  2 
v  =  n  +  V/  . 

Thus,  according  to  (1),  when  an  alert  occurs 


(23) 


V 


2 


CR-R0)2 


(24) 


Then  substituting  from  (24)  into  (15)  leads  to  the  result 


R(R-R  )  -  i/d2(R-R  )2  +  (D2-R2)VflT‘ 

Of  O  _  _ O 


(R-R0)2  +  V02t2 


T  . 


(25) 


The  following  argument  shows  that  when  the  aircraft  are 
closing  tp  is  an  increasing  function  of  |V0|.  Define  v  by 


(26) 


Then  (25)  can  be  written 


|V0|t  =  (R-R0)v. 
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For  closing  aircraft,  when  there  is  an  alert,  R  >  RQ  according 
to  (1).  Also,  if  the  separation  standard  has  not  already  been 
violated  R  >  D.  Then  it  follows  by  inspection  of  (27)  that  t^ 
is  an  increasing  function  of  v  and  hence,  by  (26),  an  increas¬ 
ing  function  of  |VQ|  as  stated. 

This  fact,  combined  with  (18),  (19)  and  (24),  provides  the 
range  of  VQ  values  for  which  a  maneuver  is  necessary,  after  an 
alert,  in  order  to  maintain  the  desired  separation.  That  is, 
the  miss  distance  will  be  at  least  D  if  and  only  if  Vq  lies 
outside  the  range  defined  by 


(R-R  )D 

os  • 


(28) 


The  probability  that  VQ  is  in  the  range  defined  by  (28)  is 
given  by 

2 


P ( 0  s  |Vj  £  Vn)  = 


0  *  2°2 
e  dV. 


(29) 


The  probability  that  the  alert,  having  occurred,  is  un¬ 
necessary  is  the  complement: 


pu  -  1  -  s  IV  s  V 


(30) 


The  probability  Py  as  a  function  of  the  separation  R  that 
exists  when  an  alert  occurs  is  illustrated  in  Fig.  5  for  various 
values  of  the  DMOD  parameter  R  .  For  these  cases  the  miss  dis¬ 
tance  standard  D  is  1,000  ft  and  the  quantity  ax  is  4,325  ft, 
which  corresponds  to  a  a  of  173  ft/sec  and  a  x  of  25  sec. 

Figure  5  indicates  that  for  aircraft  separations  beyond 
0.5  nmi  or  so,  reducing  DMOD  has  a  relatively  small  effect  on 
reducing  the  probability  that  an  alert  is  unnecessary  and  hence 
on  reducing  the  number  of  unnecessary  alerts.  On  the  other  hand. 
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FIGURE  5.  Probability  that  an  alert,  having  occurred 
unnecessary. 


reducing  DMOD  has  a  substantial  effect  on  the  reduction  of  un¬ 
necessary  alerts  only  for  separations  less  than  1,500  ft  or  so. 

This  phenomenon  has  an  interesting  parallel,  according  to 
the  remark  in  Section  IV-B  concerning  the  change  ARQ  in  DMOD 

that  is  equivalent  to  a  change  At  in  t.  The  point  was  made 
AR  AR 

that  At  is  -  — ■— ,  which  is  equal  to  -b—  t  when  an  alert  occurs 
R  R-K0 

implying  that  the  equivalent  change  in  t  is  largest  when  the 
separation  R  is  near  the  value  RQ  of  DMOD. 

Given  that  an  alert  has  occurred,  a  measure  of  the  algo¬ 
rithm  safety  is  the  probability  distribution  P(0  £  tQ  £  t)  of 
the  available  maneuvering  time  t^.  If  VQ  in  (29)  is  replaced 
by  the  more  general  | VQ | ,  the  probability  distribution  for  tD 
can  be  obtained  in  a  straightforward  way  from  (29)  with  the  aid 
of  (26)  and  (27).  Curves  illustrating  this  probability  distri¬ 
bution  are  presented  in  Fig.  6  for  various  values  of  the  DMOD 
parameter  Rq  ft  and  values  of  the  separation  R  ft  for  which 
DMOD  has  a  significant  effect  on  the  number  of  unnecessary 
alerts,  again  with  assumed  values  of  1,000  ft  for  the  miss  dis¬ 
tance  standard  D  and  A, 325  ft  for  ax. 

It  can  be  seen  that  each  probability  distribution  curve  in 
Fig.  6  starts  from  0  at  sec,  given  by  (20),  rises  to  a 


value  1-Py  given  by  (29),  at  tDmax  sec,  given  by  (21),  and  then 
becomes  horizontal.  Between  tDmin  and  tDmax  the  curves  are 
convex;  hence  a  lower  bound  on  the  probability  distribution  for 
tp  can  be  obtained  by  replacing  the  curve  in  this  interval  by  a 
straight  line,  which  is  indicated  by  a  dashed  line  in  Fig.  6. 
Clearly,  this  lower  bound  can  be  constructed  by  inspection  of 
the  curves  in  Figs.  4  and  5  with  no  need  for  further  calcula¬ 
tion  other  than  the  obvious  changes  in  scale  that  are  required 
for  consistency  with  the  measurement  units  adopted  in  Fig.  6. 


*The  implication  is  that  any  additional  alerts  suppressed  by 
the  reduction  of  DMOD  must  have  been  necessary,  i.e.,  their 
suppression  will  result  in  a  violation  of  the  miss  distance 
standard. 
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FIGURE  6.  Probability  distribution  for  the  time  to  a  separation 
equal  to  the  miss  distance  standard  when  an  alert 
occurs . 


E.  DISCUSSION  OF  SAFETY  VERSUS  UNNECESSARY  ALERTS 


An  examination  of  Pigs.  5  and  6  reveals  some  paradoxical 
relationships  between  the  aircraft  separation  R  at  the  time  of 
an  alert,  the  value  of  the  effective  DMOD  parameter  Ro>  the 
miss  distance  standard  D,  the  probability  that  the  alert  is  un¬ 
necessary,  and  the  probable  amount  of  time  available  for  a  maneu¬ 
ver  when  the  alert  is  necessary.  These  relationships  not  only 
pertain  to  the  utility  of  the  collision  algorithm  in  providing 
adequate  maneuver  time,  but  they  may  also  affect  the  aircraft 
pilot's  subjective  evaluation  of  the  reliability  of  the  BCAS. 

According  to  Fig.  5,  for  example,  when  the  effective  DMOD 
is  greater  than  or  equal  to  the  miss  distance  standard  the 
probability  that  an  alert,  having  occurred,  is  unnecessary  is 
always  high.  Moreover,  the  smaller  the  aircraft  separation  at 
the  time  of  the  alert,  the  more  likely  it  is  that  the  alert  is 
unnecessary.  To  some  extent  this  is  also  the  case  for  DMOD 
values  that  are  less  than  the  miss  distance  standard  but  near 
it. 

In  addition.  Fig.  4  indicates  that,  when  an  alert  is 
necessary,  for  DMOD  values  below  but  near  the  miss  distance 
standard,  the  maneuver  time  provided  by  the  algorithm  may  be 
greater  for  smaller  separations  at  the  time  of  the  alert  than 
for  larger.  In  such  a  case,  logically,  the  encounter  should  be 
judged  more  dangerous  when  the  aircraft  separation  Is  large 
than  when  it  is  small,  a  fact  that  appears  to  violate  intuition 
and  therefore  may  be  confusing  to  the  pilot. 

Figure  6  Is  consistent  with  these  observations.  However, 
it  also  indicates  that  for  small  values  of  the  effective  DMOD 
the  collision  algorithm  is  not  only  almost  worthless  unless  the 
separation  between  the  conflicting  aircraft  is  large,  but  may 
actually  introduce  additional  subjective  dangers.  For  example, 
the  curve  labeled  R  -  0  ft,  R  ■  1,100  ft  shows  that,  for  the 
assumed  parameter  values,  the  a  priori  probability  that  an  alert 
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having  occurred,  is  unnecessary  is  nearly  60  percent.  This  is 
for  a  separation  that  is  within  10  percent  of  the  allowable  miss 
distance.  On  the  other  hand,  if  the  alert  should  happen  to  be 
necessary  the  same  curve  indicates  that  the  maneuver  time  made 
available  by  the  collision  algorithm  will  never  be  more  than 
4.5  sec  and  the  probability  is  quite  high  that  it  will  be  less 
than  3  sec.*  Therefore,  in  this  case  the  collision  algorithm 
is  totally  counterproductive:  it  produces  a  lot  of  false 
alerts  and  never  provides  adequate  maneuver  time  when  it  gener¬ 
ates  a  true  alert.  For  DMOD  values  greater  than  zero,  but 
small  compared  with  the  miss  distance  standard,  the  same  phenom¬ 
enon  can  be  observed. 

Figure  7,  which  gives  the  probability  that  an  alert  is 
necessary  as  a  function  of  r  for  a  miss  distance  standard  of 
1,000  ft,  indicates  that  the  unnecessary  alert  problem  is 
relieved  rapidly  with  decreasing  t  for  t  values  smaller  than 
25  sec.  However,  Fig.  4  and  Fig.  6  taken  together  show  that 
the  maneuver  time  provided  by  the  algorithm  decreases  in  the 
same  proportion  as  t. 

F.  NON  CO-ALTITUDE  ENCOUNTERS 

1. 

When  two  encountering  aircraft  are  at  different  altitudes 
their  measured  separation  is  actually  a  slant  range  Rs  given  by 


It  should  be  remembered  that  the  probability  distributions 
shown  in  Fig.  6  are  a  priori.  Once  it  is  granted  that  an  alert 
is  necessary,  the  probability  scale  changes  so  that,  in  fact, 
the  horizontal  line  for  each  curve  is  positioned  at  a  prob¬ 
ability  level  equal  to  one. 


The  Effect  of  Using  Slant  Range  Instead  of  Horizontal  Range 
in  the  Horizontal  Collision  Algorithm 


where  z  is  the  relative  altitude  of  one  with  respect  to  the 
other  and  R  is  the  horizontal  separation.  On  the  other  hand, 
the  collision  algorithm  (1)  is  a  hazard  criterion  that  is  based 

on  the  horizontal  separation  and  separation  rate  of  the  air- 

•  • 

craft.  Thus,  using  R<,  in  place  of  R  (and  R^  in  place  of  R)  is 
tantamount  to  incorporating  additional  measurement  errors  in 
the  criterion. 

The  effect  of  this  additional  error  is  compensated  to  some 
extent  by  the  fact  that  there  already  exists  some  vertical 
separation  between  the  encountering  aircraft,  and,  therefore, 
less  maneuver  time  is  needed  to  achieve  the  specified  vertical 
separation  after  the  BCAS  generates  a  positive  command.  However, 
there  is  no  comparable  effect  to  mitigate  the  influence  of  the 
additional  error  on  increasing  the  chance  that  the  collision 
algorithm  will  fail  to  provide  a  necessary  alert.  Thus,  it 
seems  useful  to  examine  some  of  the  implications  of  ignoring  the 
difference  between  R<,  and  R  when  encountering  aircraft  are  not 
co-altitude. 

First  of  all,  it  should  be  recognized  that  the  effective 
error  induced  in  the  separation  rate  R  by  treating  the  slant 
range  as  if  it  were  the  horizontal  separation  has  a  special 
functional  dependence  on  R  and  R.  From  (30)  It  follows  that, 
for  a  constant  relative  altitude  z. 


Thus,  the  effective  error  in  the  separation  and  in  the 
separation  rate  due  to  using  slant  range  in  place  of  the  hori¬ 
zontal  separation  are  given  by 


When  the  aircraft  are  closing,  R  is  negative;  hence, 
according  to  (32),  eps  and  eps  are  both  negative.  If  the 
measured  slant  range  R<,  is  used  in  (1)  in  place  of  R,  an  alert 
occurs  when 

R  “  ERS  +  T  ^R  eRS  ^  <  Ro* 

i.e.,  when 

R  +  tR  <  R0  +  aRS  +  t eps.  (33 

The  implication  of  (33)  is  that  use  of  the  slant  range  in 
place  of  the  horizontal  separation  in  effect  reduces  the  DMOD 
parameter  by  the  magnitude  of  the  total  effective  error  e<, 
defined  by 

eS  "  eRS  +  teRS* 

The  error  es  is  inherently  more  serious  than  real  measurement 
errors  because  it  always  reduces  the  protection  provided  by 
DMOD,  whereas  measurement  errors  may  actually  increase  the 
protection  on  occasion. 

The  hazard  region  diagram  of  J.  Holt,  an  example  of  which 
for  co-altitude  aircraft  is  illustrated  in  Fig.  1,  can  furnish 
some  insight  concerning  the  sensitivity  of  the  protection  pro¬ 
vided  by  the  collision  algorithm  to  altitude  differences  be¬ 
tween  encountering  aircraft.  The  diagram  in  the  case  of  non 

co-altitude  aircraft  represents  the  region  defined  by  (1)  when 

•  • 

Rs  is  substituted  for  R  and  R<,  for  R. 

The  result  of  this  substitution  and  further  substitutions 
from  (2)  and  (31)  is  the  equation 

Rs  =  cos  9  +  RQ  (35! 

for  the  hazard  region  boundary.  The  curve  is  given  by  (35)  in 
the  same  polar  coordinate  system  used  previously  for  the  case 


of  co-altitude  aircraft.  The  quantity  Rs  is  given  by  (30), 
with  R  interpreted  as  the  radial  coordinate  and  z  as  a  fixed 
parameter.  In  (35)  the  quantities  t,  V  and  Rq  are  also  fixed 
parameters,  and  0  is  the  angular  coordinate. 

For  each  of  the  limiting  cases  V  =  0  and  Rq  =  0  the  hazard 
region  boundary  reduces  to  a  circle.  For  V  =  0  the  equation 
for  the  boundary  is 

R2  =  Rq2  -  z2;  (36) 

for  R  a  0  it  is 
o 

R2  -  tVR  cos  0  +  z2  =  0.  (37) 

According  to  (36),  the  circle  corresponding  to  V  *  0  is 
centered  at  the  origin  of  the  coordinate  system  and  has  a  radius 

V2  2 

R  -  z  .  According  to  (37),  the  circle  correspond- 

°  tV 

ing  to  Rq  ■  0  is  centered  on  the  polar  axis  a  distance  —  to 

the  right  of  the  origin  of  the  coordinate  system  and  has  a 

VtV  P  P 
(-y-)  -  z  . 

In  either  case  all  protection  disappears  if  the  circle 
radius  is  zero.  Thus,  when  the  relative  velocity  is  zero  there 
will  be  no  protection  if  the  difference  in  aircraft  altitudes 
is  greater  than  or  equal  to  DMOD.#  This  is  also  true  for  the 
case  of  zero  DMOD#  if  the  aircraft  altitude  difference  is 

tV 

greater  than  or  equal  to  the  quantity  —  . 

Thus,  in  the  previous  example  wherein  the  DMOD  was  assumed 
to  be  824  ft  and  x  was  assumed  to  be  25  sec,  all  protection 
vanishes  when  the  altitude  difference  is  no  smaller  than  824  ft. 
If  DMOD  is  zero  and  the  relative  velocity  is  50  ft/sec,  the 

#In  this  context  the  reference  is  to  the  effective  DMOD,  which 
is  the  remainder  after  subtracting  the  total  measurement  error 


protection  would  not  exist  for  aircraft  whose  altitude  differ¬ 
ence  is  as  large  as  625  ft. 


Figures  8,  9>  10  and  11  illustrate,  in  terms  of  hazard 
region  diagrams,  the  effect  of  altitude  differences  on  protec¬ 
tion  by  the  collision  algorithm  for  various  values  of  Vt  and 
altitude  difference  z  measured  in  units  of  Rq  ft.  As  might  be 
expected  from  consideration  of  the  two  extreme  cases,  V  ■  0 
and  Rq  *  0,  much  protection  is  lost  when  the  altitude  difference 
is  a  significant  fraction  of  DMOD  and  the  product  Vx  is  small 
compared  to  R  . 


Projected  Vertical  Miss  Distance 


In  the  case  of  non  co-altitude  aircraft  encounters  the 
active  BCAS  logic  (Ref.  5)  requires  that  the  relative  altitude 
rate  of  the  encountering  aircraft  be  projected  ahead  so  as  to 
obtain  an  estimate  of  the  vertical  miss  distance  at  the  time 
of  minimum  horizontal  separation.  The  time  used  for  this  pur- 

•o 

pose  is  the  ratio  -  -r  (referred  to  as  "true  tau" )  unless  the 

R 

ratio  exceeds  a  pre-selected  constant,  called  TVPCMD,  in  which 
event  it  is  that  constant  that  is  used. 


It  would,  perhaps,  be  more  useful  to  obtain  the  vertical 
miss  distance  that  occurs  at  the  time  when  it  is  estimated  that 
the  horizontal  separation  standard  will  be  violated.  This  is 
the  time  t^  given  by  (15)-  Since  t^  depends  upon  the  magnitude 

V  of  the  relative  velocity  and  that  quantity  is  not  known  in 
general,  it  would  be  necessary  to  use,  instead,  a  value  taken 
from  the  time  interval  between  ,  given  by  (20),  and  tDmax> 

given  by  (21),  however.  The  safest  value  for  this  purpose  is 

tDmax*  of  course- 


The  use  of  t^max  in  place  of  true  tau,  instead  of  the  con¬ 
stant  TVPCMD,  when  the  horizontal  separation  rate  is  small  also 


ensures  a  low  probability  of  missing  a  threat.  The  alert  rate 


will  tend  to  be  even  smaller  because  tDmax 
true  tau.  This  can  be  seen  as  follows. 


is 


always 


less  than 
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Z  AND  HORIZONTAL  COORDINATES  MEASURED  IN  UNITS  OF  R0  FT 


FIGURE  10.  The  hazard  region  defined  by  a  collision  algorithm 

when  slant  range  is  used  instead  of  horizontal  range 


Z  AND  HORIZONTAL  COORDINATES  MEASURED  IN  UNITS  OF  R0  FT 


According  to  (1)  at  the  time  of  a  horizontal  alert 


R  Rt  . 
R  “  R'Ro 


(38) 


On  the  other  hand,  according  to  (20)  and  (21) 


f  -  R2-d2  _ 

cDmax  "  R(R-RQ) 


Now,  it  is  certainly  true  that 


(39) 


(40) 


Then,  upon  multiplying  (40)  by  t  and  dividing  by  R-RQ  it  will 
be  seen  to  follow  from  (38)  and  (39)  that 


4-  <  —  ii 

rDmax  £  » 

as  asserted. 

G.  EFFECTIVE  HAZARD  REGION  DEFINITION 


(41) 


Reference  12  contains  the  logical  basis  for  selecting  col¬ 
lision  algorithm  parameters  when  the  only  available  data  is 
that  assumed  for  active  BCAS,  namely,  vertical  and  horizontal 
separations  and  separation  rates.  There  it  is  demonstrated 
for  horizontal  criteria  that  the  modified  tau  algorithm  (1)  is 
the  most  efficient  one  possible  in  the  sense  that  it  is  both 
necessary  and  sufficient  to  protect  against  all  encounters  for 
which  a  priori  bounds  on  relative  acceleration  and  measurement 
errors  and  a  required  time  to  perform  an  escape  maneuver  are 
assumed . 

When  the  parameters  R  and  x  appearing  in  (1)  are  chosen 
properly  there  is  no  need  for  an  additional  minimum  range 
criterion.*  Conversely,  the  exclusion  of  any  part  of  the 


hazard  region  defined  by  (1)  introduces  the  risk  that  a  neces¬ 
sary  alert  will  be  missed. 


The  same  argument  that  supports  the  use  of  the  horizontal 
algorithm  (1)  in  connection  with  horizontal  separation  data 
also  applies  to  a  similar  algorithm,  such  as  that  employed  in 
active  BCAS  logic,  in  connection  with  vertical  separation  data. 
Therefore,  a  vertical  collision  algorithm  should  have  the  same 


form  as  (1),  with  a  DMOD  parameter  analogous  to  Rq  to  account 


for  measurement  error  and  vertical  acceleration.  Because  of 


its  meaning  the  value  of  t  should  be  the  same  as  that  used  in 


the  horizontal  alert  criterion. 


In  the  present  version  of  active  BCAS  logic  (Ref.  5)  both 
the  horizontal  and  vertical  collision  algorithms  deviate  in 
certain  respects  from  these  precepts.  In  view  of  the  principles 
established  with  mathematical  precision  by  Ref.  12,  it  is  dif¬ 
ficult  to  understand  how  these  deviations  can  be  justified. 


The  active  BCAS  vertical  algorithm  has  the  form 


z  +  rz  <  0,  z  <  0 


(42) 


z  <  z  .  z  >  0 
o 


where  z  is  the  vertical  separation,  t  is  the  vertical  separa¬ 
tion  rate,  and  zQ  is  a  fixed  minimum  vertical  separation  thresh¬ 
old.  Using  (42)  instead  of  a  modified  tau  criterion  is  equiva¬ 
lent  to  using,  instead  of  (1),  a  horizontal  algorithm  with  zero 
DMOD  along  with  a  minimum  range  criterion.  This  will  not  only 
result  In  missed  necessary  alerts  on  some  occasions  but  will 
also  generate  unnecessary  alerts  on  others.* 


Actually,  the  active  BCAS  logic  does  provide  for  the  removal 
of  some,  but  not  all,  of  the  unnecessary  hazard  space  when 
there  is  a  positive  vertical  separation  rate. 


a 


In  the  case  of  the  horizontal  algorithm,  the  active  BCAS 
logic  removes  part  of  the  hazard  space  when  there  is  a  positive 
horizontal  separation  rate.*  According  to  Ref.  5  this  is  done 
in  order  to  reduce  the  amount  of  time  that  an  alert  condition 
remains  in  effect  when  the  encountering  aircraft  are  separating 
The  implication  is  that  the  risk  is  less,  or  becomes  less  in 
due  course,  during  an  alert  that  occurs  when  encountering  air¬ 
craft  are,  or  appear  to  be,  separating  rather  than  closing. 

However,  it  is  possible  to  dispute  this  assumption  on  two 
accounts.  The  algorithm  (1)  provides  a  hazard  region  even  when 
there  is  a  positive  separation  rate  R  because  of  the  presence 
of  the  DMOD  parameter  Rq,  which  was  introduced  in  order  to  com¬ 
pensate  for  measurement  error  and  some  relative  acceleration. 
Therefore,  the  assumption  may  be  incorrect  because  (a)  due  to 
measurement  error,  the  aircraft  may  actually  be  closing  even 
though  they  appear  to  be  separating,  and  (b)  if  the  aircraft 
are  initially  separating  slowly  while  traveling  on  nearly  par¬ 
allel  flight  paths  a  turn  by  one  of  them  after  the  alert  should 
have  occurred  could  rapidly  change  the  separation  rate  from 
positive  to  negative. 

In  fact,  neglecting  an  alert,  or  prematurely  eliminating 
an  alert  condition,  may  be  more  hazardous  In  the  region  of  posi¬ 
tive  separation  rates  than  elsewhere  because  in  that  region  the 
separation  R  is  necessarily  small.  This  will  be  the  case  if 
the  effective  DXCD  parameter  Rq  is  less  than  the  miss  distance 
standard  D  and  the  separation  R  is  near  D,  as  Figs.  4  and  5 
indicate . 


Chapter  II,  Section  B. 


V.  SUMMARY  OF  CONCLUSIONS 


A.  VALIDITY  OF  A  MATHEMATICAL  MODEL  OF  AIR  TRAFFIC 

A  mathematical  model  for  low-altitude  (below  10,000  ft) 
air  traffic,  originally  formulated  to  fit  the  Los  Angeles  Stand¬ 
ard  Traffic  Model,  has  been  tested  against  empirically-derived 
traffic  statistics  for  the  Houston  Intercontinental  air  terminal 
region.  According  to  the  theory,  the  horizontal  cartesian 
coordinates  of  an  aircraft,  relative  to  some  geographical  center, 
are  identically  distributed,  independent  Gaussian  random  vari¬ 
ables  with  mean  zero  and  a  common  standard  deviation;  the  same 
statistical  assumptions  are  postulated  for  the  aircraft's  veloc¬ 
ity  components. 

The  theoretical  model  is  consistent  with  processed  empiri¬ 
cal  data  derived  by  the  Mitre  Corporation  from  ARTS  tapes 
containing  radar  records  of  65  hours  of  air  traffic  around  the 
Houston  terminal.  Mitre  processed  the  data  to  expose  and  study 
mid-air  conflicts  that  would  be  defined  by  active  BCAS  alert 
logic  at  various  parameter  settings.  Results  that  are  compared 
here  with  the  theory  are  of  two  types.  One  is  in  the  form  of 
histograms  showing,  for  aircraft  in  conflict:  (1)  the  dis¬ 
tribution  of  absolute  speeds,  (2)  the  distribution  of  closing 
speeds,  (3)  the  distribution  of  separations.  The  other  is 
the  reduction  that  occurs  in  the  number  of  alerts  when  BCAS 
logic  parameters  are  reduced  by  certain  specified  amounts. 

The  theoretical  traffic  model  fits  the  histogram  data  best 
when  the  velocity  component  standard  deviation  is  taken  to  be 
106  kn.  The  accuracy  of  the  fit  appears  to  be  independent  of 
the  model's  traffic  density  standard  deviation. 


The  theory  also  predicts  the  alert  reductions  that  were 
observed  by  Mitre  when  BCAS  logic  parameters  were  changed.  The 
predictions  are  virtually  independent  of  the  models’  traffic 
density  standard  deviation  and  are  accurate  for  three  different 
values  of  its  velocity  component  standard  deviation  that  were 
tried:  72.2  kn  (the  Los  Angeles  Standard  Traffic  Model  value), 

83*3  kn  (the  99  percent  confidence  level  that  no  aircraft 
violates  the  250  kn  speed  limit),  and  106  kn  (the  value  that 
provides  the  best  fit  to  the  Mitre  histogram  data).  This  leads 
to  the  conjecture  that  alert  rate  changes  due  to  changes  in 
BCAS  logic  parameters  are  essentially  independent  of  the  air 
traffic  environment. 

B.  BCAS  ALERT  LOGIC  PROTECTION 

1 .  Alert  Logic  Parameter  Selection 

a.  Some  allowance  for  acceleration  should  be  made  in 
selecting  the  DMOD  parameter  because,  aside  from  the  possibility 
of  turns,  certain  aircraft  may  not  be  able  to  perform  vertical 
avoidance  maneuvers  without  introducing  an  additional  horizon¬ 
tal  velocity  component  into  the  forward  motion. 

b.  A  DMOD  parameter  should  be  included  in  the  vertical 
tau  alert  criterion  because  of  altitude  measurement  errors  and 
vertical  acceleration. 

c.  In  desensitization  a  reduction  of  x  is  preferable,  in 
principle,  to  a  reduction  in  DMOD  because  reducing  DMOD  reduces 
protection  for  all  aircraft,  whereas  reducing  x  may  not  affect 
aircraft  that  can  react  quickly  to  an  alert. 

d.  If  reasonable  x  values  (see  Pig.  3)  are  maintained  when 
alert  logic  parameters  are  reduced,  the  probability  of  missing 
an  alert  that  is  necessary  to  avoid  violating  an  adequate  mini¬ 
mum  miss  distance  standard  will  be  low. 

e.  The  projected  vertical  miss  distances  rule  presently 
used  in  active  BCAS  logic  relies  upon  the  substitution  of  a 


pre-selected  constant,  TVPCMD,  for  true  tau,  which  is  the  sepa¬ 
ration  divided  by  the  closing  rate,  when  the  value  of  true  tau 
exceeds  that  constant.  A  better  way  to  deal  with  large  values 
of  true  tau  (i.e.,  small  separation  rates)  might  be  to  use 
tDmax*  Siven  by  (21),  in  place  of  true  tau. 

2.  Protection  And  The  Alert  Rate 


a.  For  a  given  minimum  miss  distance  standard  D  and  a 
separation  R  when  there  is  an  alert,  either  the  alert  is  un¬ 
necessary  or  the  amount  of  available  maneuver  time  lies  between 

tDmin»  ®iven  (2°)>  and  ^max*  ®lven  (21>* 

b.  When  the  effective  DMOD  is  greater  than  or  near  the 
minimum  miss  distance  standard,  the  smaller  the  aircraft  separa 
tion  at  the  time  of  an  alert,  the  more  likely  it  is  that  the 
alert  is  unnecessary. 

c.  When  the  effective  DMOD  is  smaller  than  but  near  the 
minimum  miss  distance  standard  the  maneuver  time  available 
after  an  alert  may  be  greater  for  smaller  separations  of  the 
aircraft  in  conflict  than  for  larger.  This  may  violate  a 
pilot’s  intuition  and  tend  to  reduce  his  acceptance  of  the 
system  unless  he  is  forewarned  of  the  paradox. 

d.  The  possibility  of  pilot  disenchantment  with  the  sys¬ 
tem  is  enhanced  by  the  fact  that  the  percentage  of  unnecessary 
alerts  will  always  be  high.  In  fact,  this  is  true  even  for 
small  values  of  DMOD,  for  which  the  available  maneuver  time 
after  a  necessary  alert  will  generally  be  too  small  to  be  use¬ 
ful.  This  will  be  true  unless  the  separation  is  large,  indicat¬ 
ing  a  head-on  or  near  head-on  collision.  Therefore,  it  may  be 
better  to  turn  the  system  off  than  to  desensitize  it  to  very 
small  DMOD  values. 

3 .  Dubious  Compromises  of  Protection 

a.  The  use  of  slant  range  in  place  of  a  horizontal  range 
estimate  may  seriously  compromise  protection  because  the  effec 
tive  error,  unlike  random  measurement  errors,  always  reduces 
protection. 


b.  When  aircraft  are  apparently  separating  instead  of 
closing,  reducing  the  hazard  region  defined  by  BCAS  alert  logic 
which  is  the  current  practice,  is  dangerous  because  of  measure¬ 
ment  error  and  possible  unforeseen  accelerations  and  the  fact 
that  in  such  a  region  the  aircraft  separation  is  necessarily 
small. 
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ANALYSIS  AND  APPLICATION  OF  A  THEORETICAL 
MODEL  OF  AIR  TRAFFIC  STATISTICS 

A.  VELOCITY  DISTRIBUTIONS 

As  observed  in  Ref.  A-l,  low  altitude  aircraft  velocities 
in  the  Los  Angeles  Standard  Traffic  Model  (Ref.  A-2)  appear  to 
be  normally  distributed  with  no  directional  bias.  More  speci¬ 
fically,  the  two  horizontal  velocity  components  associated  with 
every  aircraft  appear  to  be  independent  random  variables  with 
identical  Gaussian  probability  distributions,  each  with  zero 
mean  and  a  standard  deviation  of  122  ft/sec  (72.2  kn).  This 
implies  that  the  aircraft  speeds  are  Rayleigh  distributed,  i.e., 
have  a  probability  density  of  the  form 


where  S  is  the  aircraft  speed  and  in  this  case  the  parameter 
oQ  has  the  value  122  ft/sec  or  72.2  kn. 

It  also  follows  that  the  components  of  the  relative  veloc¬ 
ity  between  two  aircraft  selected  at  random  are  independent, 
identically  distributed  Gaussian  random  variables  with  zero 
mean  but  with  a  standard  deviation  av  of  /IT  times  122  ft/sec, 
or  173  ft/sec.  Thus,  the  aircraft  separation  rate,  because  it 
is  equal  to  the  radial  component  of  their  relative  velocity, 
is  a  Gaussian  distributed  random  variable  with  zero  mean  and  a 
standard  deviation  of  173  ft/sec. 


Reference  A- 3  contains  a  figure  with  a  histogram  depicting 
the  percentage  distribution  of  speeds  of  aircraft  in  conflict 
in  the  region  around  Houston  Intercontinental  Airport.  Actually 
the  figure  contains  two  such  histograms,  for  one  of  which  the 
definition  of  a  conflict  is  based  on  a  fixed  set  of  parameters 
in  the  BCAS  alert  logic,  while  for  the  other  the  parameters  were 
permitted  to  vary  in  accordance  with  the  local  desensitization 
zone  rules.  Only  the  data  generated  by  the  fixed  parameter 
logic  is  considered  here. 

It  is  assumed  here  that  the  number  of  aircraft  in  conflict 
above  10,000  ft  is  negligible  during  the  65-hour  period  of  ARTS 
tape  data  analyzed  by  the  Mitre  Corporation  for  the  program 
summarized  in  Ref.  A-3.  This  is  tantamount  to  an  assumption 
that  restricting  the  data  to  aircraft  in  conflict  is  equivalent 
to  taking  a  sample  of  low-altitude  aircraft,  albeit  one  that 
is  somewhat  biased  against  low  speeds. 

With  these  assumptions  the  histogram  supplied  by  Ref.  A-3 
can  be  used  to  ascertain  whether  a  probability  density  of  the 
form  ( A- 1 ) ,  deduced  originally  from  the  Los  Angeles  Standard 
Traffic  Model,  will  fit  the  Houston  traffic  environment  as  well. 
For  this  purpose  It  should  be  observed  that  the  probability 
distribution  corresponding  to  the  density  defined  by  (A-l)  is 
given  by 
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Figure  A-l  reproduces  the  aircraft  speed  histogram  of 
Ref.  A-3  superposed  on  three  probability  distribution  curves 
calculated  by  means  of  (A-2).  The  dashed-line  curve  corres¬ 
ponds  to  the  Los  Angeles  Standard  Traffic  Model  cq  value  of 
72.2  kn  (122  ft/sec).  The  solid-line  curve  corresponds  to  a 
oQ  value  of  83.3  kn,  which  would  imply  a  99  percent  confidence 
level  that  no  aircraft  violates  the  speed  limit  of  250  kn.  The 


remaining  curve  corresponds  to  a  oq  value  of  106  kn  and  is 
included  for  reasons  that  will  become  apparent  later. 

It  is  evident  from  the  figure  that  the  Los  Angeles  Traffic 
Model  distribution  does  not  fit  the  histogram  as  well  as  the 
other  two.  Although  neither  of  those  could  be  regarded  as 
incompatible  with  the  empirical  data,  it  is  clear  that  a  dis¬ 
tribution  corresponding  to  a  aQ  value  between  83-3  and  106  kn 
would  fit  best. 

B.  TRAFFIC  DENSITY  DISTRIBUTION  AND  THE  ALERT  PROBABILITY 

Reference  A-3  does  not  provide  sufficient  data  to  infer 
the  horizontal  or  vertical  distribution  of  aircraft  in  the 
Houston  terminal  region.  However,  as  in  the  case  of  velocities 
statistics  from  the  Los  Angeles  Standard  Traffic  Model  can  be 
assumed. 

According  to  Ref.  A-l  low-altitude  aircraft  in  the  model 
are  distributed  about  a  traffic  center  in  such  a  manner  that 
the  horizontal  coordinates  relative  to  the  center  are  indepen¬ 
dent,  identically  distributed  Gaussian  random  variables  with 
mean  zero  and  a  standard  deviation  aR.*  The  value  of  aR  in  the 
Los  Angeles  model  is  20  nmi . 

It  follows  that  the  traffic  is  Rayleigh  distributed  in 
range  from  the  center  and  uniformly  in  direction.  It  also 
follows  that  the  separation  R  between  two  aircraft  chosen  at 
random  is  Rayleigh  distributed  with  the  probability  density 


where 


cR  =  /T  aR  . 

^Actually,  the  distribution  given  by  Ref.  A-l  is  bimodal,  but 
one  mode  is  so  dominant  that  the  other  can  be  neglected  with¬ 
out  affecting  any  results  obtained  here. 
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This,  together  with  the  assumed  Gaussian  distribution  for 
the  separation  rate  A,  can  be  used  to  calculate  the  horizontal 
alert  probability.  In  terms  of  t  and  the  DMOD  parameter  Rq 
the  condition  to  be  imposed  for  this  purpose  is 

R  +  t  ft  *  Rq*. 

The  probability  that  two  aircraft  selected  at  random  will 
be  in  a  horizontal  alert  condition  can  be  expressed,  in  accord¬ 
ance  with  ( A—  3 ) >  by 
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This  expression  can  be  transformed,  by  standard  means,  into  one 
involving  only  elementary  functions  and  normal  distribution 
functions  that  are  easy  to  evaluate  numerically. 


First,  however,  it  is  advantageous  to  simplify  ( A—  4 )  by 

A 

changing  variables  so  that  the  parameters  RQ,  a^,  t,  and  oR  are 
replaced  by  dimensionless  quantities.  It  is  found,  in  fact, 
that  the  introduction  of  quantities  p  and  <  defined  by 


o  r 


(A-5) 
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#The  two  other  horizontal  threat  criterion  conditions  given  on 
p.  5  should  have  a  negligible  effect  on  any  conclusions  based 
on  the  analysis  of  this  appendix. 
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reduces  the  number  of  Independent  parameters  from  four  to  two. 
The  result  is  an  expression  equivalent  to  (A-4): 


P(R0,crv,T,ap)  =  P(k,p) 
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The  function  I2(k,p)  in  (A-7)  can  be  simplified  further  by 
completing  the  square  of  the  exponent  in  the  integrand;  i.e.. 
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Figure  A-2  provides  curves  depicting  the  probability  P(k,p) 
as  a  function  of  p  for  various  values  of  k.  The  range  of  <  and 
p  in  the  figure  should  be  sufficient  to  cover  all  combinations 

A 

of  values  of  the  original  parameters,  R  ,  ay,  x,  aR  that  are 
likely  to  be  of  interest  for  analyses  of  BCAS.  The  adequacy  of 
the  range  of  <  and  p  can  be  verified  by  substituting  numerically 
typical  values  for  R  ,  av>  x,  and  aR  into  (A-5);  e.g., 

■  1  nmi, 
o 

cy  ■  100  kn. 
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FIGURE  A- 2 


Probability  that  two  aircraft  selected  at  random 
will  satisfy  a  horizontal  alert  condition. 
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According  to  (A-5),  if  the  DMOD  parameter  Rq  is  changed  by 
a  factor  a  the  quantity  p  is  changed  by  the  same  factor  but  ic 
is  unaffected.  However,  if  t  is  changed  by  a  factor  a  the  quan¬ 
tity  <  is  changed  by  the  same  factor  and  p  by  the  factor  1/a; 
in  fact,  the  product  tcp  is  independent  of  t  and  av- 

C.  OBSERVED  EFFECT  OF  PARAMETER  CHANGES  IN  THE  ALERT  LOGIC 

Reference  A-3  gives  estimates  of  the  percentage  change  in 
the  number  of  alerts  for  the  Houston  traffic  data  due  to  per¬ 
centage  changes  in  t  and  in  DMOD  (RQ).  One  would  expect  that 
the  probability  depicted  by  the  curves  in  Fig.  A-2  will  be  pro¬ 
portional  to  the  number  of  alerts  over  any  given  time  period 
and  therefore  that  the  same  percentage  changes  should  be  observed 
in  that  probability.  If  so,  it  might  be  expected  that  the 
results  reported  in  Ref.  A-3,  when  compared  with  similar  changes 
induced  in  the  probability  depicted  in  Fig.  A-2,  would  provide  a 
means  for  estimating  the  quantities  and  av  for  the  Houston 
traffic. 

Reference  A-3  reports: 

(1)  decreasing  t  from  25  to  20  sec  reduces  the  number  of 
alerts  by  31-5  percent; 

(2)  decreasing  DMOD  (Rq)  from  0.3  nmi  to  0.1  nmi  reduces 
the  number  of  alerts  by  25  percent.  However,  in  another  section 

A-10 


Ref.  A-3  states  that  for  a  370  ft  altitude  separation  standard 
(designated  as  ALIM)  for  Rq  *  0.3  nmi  and  t  =  25  sec  there  are 
181  alerts,  and  for  the  same  altitude  parameter,  when  Rq  is 
reduced  to  0.1  nmi  and  t  to  20  sec,  the  number  of  alerts  is 
reduced  to  83. 

The  partial  percentage  reductions  in  alerts  reported  by 
Ref.  A-3  are  incompatible  with  the  reported  total  reduction 
from  l8l  to  83,  which  is  equivalent  to  a  total  reduction  of 
54.1  percent.  The  combined  reported  partial  reductions  imply 
a  total  reduction  of  [1  -  (1-0. 25) (1-0. 315) ]  x  100,  or  48.6 
percent.  This  discrepancy  may  be  due  to  a  typographical  error, 
e.g.,  if  25  percent  were  inadvertently  written  in  place  of  33*1 
percent . 

A  similar  discrepancy  occurs  for  the  other  total  reduction 
figures  given  by  Ref.  A-3  for  the  case  of  a  340  ft  altitude 
separation  standard,  namely,  a  reduction  from  173  to  8l.  In 
this  case  the  total  reduction  is  53*2  percent  rather  than  the 
implied  percentage  reduction  of  48.6  percent. 

Because  of  these  discrepancies  and  the  fact  that  the  data 
concerning  the  effect  of  reducing  the  parameters  separately  is 
not  given  explicitly  in  terms  of  actual  alert  numbers  by 
Ref.  A-3,  those  reported  reductions  will  be  ignored  here.  The 
explicit  data  concerning  the  combined  effect  of  reducing  both 
t  and  DM0 D  will  be  used,  however,  and  are  displayed  in  Table  A-2. 

D.  APPLICATION  OF  THE  THEORETICAL  MODEL  TO  THE  HOUSTON  DATA 

Reference  A-3  also  provides  some  additional  data  that  can 
be  used  to  fit  the  theoretical  air  traffic  model  embodied  in 
(A-6)  to  the  observed  statistics  for  Houston.  This  data  is 
another  histogram,  reproduced  in  Fig.  A-3>  that  depicts  the 
distribution  of  aircraft  separations  for  the  conflicts  observed 
during  the  65  hours  covered  by  the  Houston  ARTS  tapes. 
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Finally,  L.  Zarrelli  has  kindly  provided  another  histogram 
that  depicts  the  distribution  of  closing  speeds  for  the  conflicts, 
It  is  reproduced  here  in  Fig.  A-4. 


For  the  data  contained  in  both  histograms  all  conflicts 
were  defined  by  level  4  alert  parameter  values:  a  t  value  of 


25  sec  and  a  DMO D  (R  )  value  of  0.3  nmi.  Similar  data,  for 


which  the  conflicts  were  determined  by  the  application  of 
desensitization  rules  according  to  the  previously  defined  local 
zones  in  which  they  occurred,  were  also  provided  by  Ref.  A-3 
and  L.  Zarrelli.  However,  because  of  the  mixing  of  parameter 
values  these  data  are  not  suited  to  the  purposes  of  this  report. 


The  theoretical  model  of  air  traffic  statistics  can  be 
tested  by  using  it  to  predict  the  percentage  changes  that  took 
place  in  the  number  of  alerts  when  certain  alert  logic  parameter 
changes  were  made  in  the  data  processing  experiment  reported  by 
Ref.  A-3.  For  this  purpose  it  should  be  remembered  that  the 


quantity  appearing  in  (A-5)  is  /Y  times  the  oq  of  the 


Rayleigh  distributions.  Thus,  the  72.2  kn  value  of  oq  becomes 


102.1  kn  and  the  83.3  kn  value  becomes  117.8  kn. 


The  quantity  tc  in  (A-5)  requires  a  knowledge  of  aR,  which 
is  associated  with  the  density  distribution  of  aircraft.  How¬ 
ever,  after  some  numerical  experimentation  it  can  be  seen  that 


the  predicted  reductions  in  alerts  are  virtually  independent 

A  « 

of  ap  Therefore,  an  arbitrarily  chosen  value  of  20  nmi  for 

/v  n* 

aR  seems  as  good  as  any  other  within  physically  reasonable 
limits . 


The  paired  values  of  p  and  ic  of  interest,  calculated  by 
means  of  (A-5),  are  those  appearing  in  Table  A-l  along  with  the 
corresponding  probability  calculations  given  by  (A-6).  Table 
A- 2  gives  the  corresponding  alert  changes  observed  empirically, 
according  to  Ref.  A-3-  Table  A-3  is  a  comparison  of  predicted 
with  observed  results. 


»No  attempt  has  been  made  in  the  present  investigation  to  explain 
this  remarkable  fact,  which  is  not  an  obvious  consequence  of 
the  general  theory. 
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THEORETICAL  DISTRIBUTION 
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Observed  percent  of  conflicts  with  closing  speeds 
less  than  X. 


TABLE  A- 1 .  PREDICTED  CHANGES  IN  THE  ALERT  PROBABILITY  BECAUSE 
OF  CHANGES  IN  ALERT  LOGIC  PARAMETERS 


x(sec) 

RQ(nm1) 

cQ(Kn) 

P  K 

Prob.  x  10^ 

%  Change 

25 

0.3 

72.2 

0.423  0.0355 

5.90 

20 

0.1 

72.2 

0.176  0.0284 

2.65 

55.1 

18 

0.1 

72.2 

0.196  0.0256 

2.22 

62.4 

25 

0.3 

83.3 

0.367  0.0409 

7.24 

- 

20 

0.1 

83.3 

0.153  0.0327 

3.89 

53.2 

18 

0.1 

83.3 

0.170  0.0294 

2.81 

61.2 

25 

0.3 

106 

0.288  0.0520 

10.4 

- 

20 

0.1 

106 

0.120  0.0416 

5.21 

49.9 

18 

0.1 

106 

0.133  0.0374 

4.30 

58.7 

TABLE  A- 

-2.  DISTRIBUTION  OF  ALERTS 

(PARTIAL  DATA) 

DMOD ( Rq ) 

ALIM 

TAU  Positive  Commands 

0.3 

nmi 

370  ft 

25  sec 

181 

0.3 

nmi 

340  ft 

25  sec 

173 

0.1 

nmi 

370  ft 

20  sec 

83 

0.1 

nmi 

340  ft 

20  sec 

81 

0.1 

nmi 

370  ft 

18  sec 

73 

0.1 

nmi 

340  ft 

18  sec 

69 

TABLE 

A-  3. 

ALERT  REDUCTION  DUE  TO  CHANGES  IN  t  AND 

DMOD 

Change  1 n 
t  From  25 
secs  to 
(sec) 

Alim 

(ft) 

Actual  % 
Change  in 
Alerts 

Predicted  % 

Change  in  Alerts  Error  % 

72.2  83.3  106  72.2  83.3 

106  °o(kn) 

20 

370 

54.1 

55.1 

53.2 

49.9 

1.0 

0.9 

4.2 

20 

340 

53.2 

55.1 

53.2 

49.9 

1.9 

0.0 

3.3 

18 

370 

59.7 

62.4 

61.2 

58.7 

2.7 

1.5 

1.0 

According  to  these  results,  for  example,  if  ■  72.2  kn 
the  predicted  reduction  in  alerts  would  be  100  x  |l  -  ^-|^jper- 
cent  or  55.1  percent.  Applied  to  the  l8l  aircraft  in  conflict 
for  an  altitude  standard  of  370  ft  this  would  yield  a  predicted 
reduction  to  81.3  aircraft,  a  result  that  differs  from  what  was 
observed  by  1.7  aircraft.  This  is  an  error  of  less  than  one 
percent . 


The  cq  value  of  83*3  kn  is  apparently  the  most  consistent 


with  the  Houston  data;  it  leads  to  the  most  accurate  prediction 
of  the  effect  of  alert  logic  parameter  changes  on  the  reduction 


of  alerts.  However,  all  of  the  oq  values  considered  lead  to 


predictions  that  are  in  excellent  agreement  with  the  observed 
results.  Thus,  it  appears  that  the  effect  of  parameter  changes 
on  the  reduction  of  alerts  is  not  very  sensitive  to  the  air 
traffic  statistics. 


The  remaining  data,  the  aircraft  closing  speed  and  separa¬ 
tion  statistics,  can  be  used  to  test  the  theoretical  traffic 


model  and  aQ  values.  For  this  purpose  it  is  necessary  to  calcu¬ 


late  the  appropriate  probabilities.  These  can  be  defined  as 
follows.  P  (x)  is  the  conditional  probability  that  two  air¬ 
craft,  known  to  be  in  conflict,  have  a  closing  speed  that  is 
less  than  x  kn,  and  PR(x)  is  the  conditional  probability  that 
their  separation  at  the  beginning  of  the  conflict  is  less  than 
x  nmi. 


Dividing  the  joint  probability  PvC(x)  that  two  aircraft 
are  in  conflict  and  have  a  closing  speed  less  than  x  by  P(k,p), 
the  probability  that  the  two  aircraft  are  in  conflict,  yields 


Pv(x) .  A  similar  calculation  using  P^^(x),  the  joint  probabil¬ 


ity  that  two  aircraft  are  in  conflict  and  have  a  separation 
less  than  x  at  the  beginning  of  the  conflict,  in  place  of  PvC(x), 
yields  PR(x) . 
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The  Joint  probabilities  can  be  written  at  once  using  the 
horizontal  collision  algorithm  condition  and  the  assumed  theo¬ 
retical  traffic  distributions.  That  is. 


PvC(x) 
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These  integrals  can  be  transformed  by  means  of  the  same 
procedure  used  in  deriving  (A-6)  from  (A-4).  The  results  are 
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where  Pc  is  given  by  (A-6).  The  required  probabilities  are 
then  given  by 
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Histograms,  constructed  by  Mitre  from  the  Houston  data, 

showing  the  observed  closing  speed  and  separation  distributions 

of  aircraft  in  conflict  are  reproduced  in  Fig.  A-4  and  A-3. 

Included  in  those  figures  are  theoretical  curves  calculated 

using  Rq  *  0.3  nmi,  r  -  25  second,  cr^  *  20  nmi#  in  (A-10)  with 

values  of  aQ  equal  to  83-3  kn  (dashed  line)  and  106  kn  (solid 

line).  In  both  cases  the  aQ  value  of  83.3  kn  produces  curves 

that  lie  above  the  corresponding  histogram.  Conversely,  the 

theoretical  distributions  based  on  a  cr  value  of  106  kn  fit 

0 

their  corresponding  histograms  reasonably  well  in  both  cases. 


Some  numerical  experimentation,  again,  demonstrates  that  the 
theoretical  distributions  are  virtually  independent  of  oR. 
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